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We present measurements of the complete spatio-temporal Fourier spectrum of Faraday waves. 
The Faraday waves are generated at the interface of two immiscible index matched liquids of different 
density. By use of a new absorption technique we are able to determine the bifurcation scenario from 
the flat surface to the patterned state for each complex spatial and temporal Fourier component 
separately. The measured surface spectra at onset are in good agreement with the predictions 
from a linear stability analysis. For the nonlinear state our measurements show in a direct manner 
how energy is transferred from lower to higher harmonics and we quantify the nonlinear coupling 
coefficients. Furthermore we find that the nonlinear coupling generates constant components in 
the spatial Fourier spectrum leading to a contribution of a non oscillating permanent sinusoidal 
deformed surface state. A comparison of a hexagonal and a rectangular pattern reveals that spatial 
resonances only can give rise to a spectrum that violates the temporal resonance conditions given 
by the weakly nonlinear theory. 

PACS numbers: 89.75.Kd, 47.54.-|-r, 47.35.-(-i, 47.20.Gv 
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1. Introduction 



The Faraday Experiment is probably the first non equi- 
librium pattern forming system that has been investi- 
gated scientifically, namely by Michael Faraday in 1831 
Nevertheless it was only recently that it was possi- 
ble to determine the complete Fourier spectrum of the 
deformed surface state 0- While an experimental anal- 
ysis of the full mode spectrum in other pattern forming 
model systems like Rayleigh-Benard or Taylor-Couette 
has been standard technique for a long time, it is the re- 
fraction of light at the free surface of a liquid that renders 
the analysis of surface waves so difficult. Quantitative in- 
formation about the patterned state up to higher orders 
is important not only to verify the validity of theoretical 
calculations 0, 3 but also to gain insight on the reso- 
nance mechanisms that form the patterns. The Faraday 
experiment is especially known for its richness of differ- 
ent patterns that are observed |E IE S IE lE^E- By 
using complex liquids [T]| , very low fill heights |l2j | or by 
introducing additional driving frequencies, highly com- 
plex ordered states like superlattices [ij, |lj] have been 
observed recently. But we will demonstrate that even 
simple patterns like lines, squares and hexagons observed 
in a single driving frequency experiment can still unveil 
unknown surprising characteristics. 



A discussion of the different attempts to reveal quan- 
titative information on the the surface elevation profile 
C{r,t) of Faraday waves is given in ^E|- The main diffi- 
culties in determining the surface elevation profile of cap- 
illary waves at the the free surface are the difference in re- 
fractive indexes of the liquid and the air, and the fact that 
the interface diffuses almost no light but rather reflects 
or transmits incoming light completely. To our knowl- 
edge there is only one optical method that overcomes this 
problem with the use Polystyrene colloids to provide light 
scatterers within the fluid [13 , but the method was only 
used in a turbulent regime. Another powerful method for 
the investigation of capillary waves on ferrofluids based 
on x-ray absorption was presented in |0 , but the related 
costs and efforts might be justified for fully opaque liq- 
uids only. To bypass the problems associated with light 
refraction and reflection on a liquid-air interface we chose 
to study the interface between two index matched liquids. 
The upper fluid is transparent, the lower one is dyed. In 
the presence of surface deformations the instantaneous 
thickness can be deduced from the intensity of the light 
transmitted trough the colored layer. From a hydrody- 
namic point of view the replacement of the air by a second 
liquid is nothing but a change of viscosity and density, 
though the low kinematic viscosity of air simplifies the 
theoretical calculations. However, the first exact theo- 
retical analysis of the linear stability problem by Kumar 
and Tuckerman ^3] was carried out for the more general 
case of a system of two layers of liquid. 
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2. Experimental setup 

The experimental setup is shown in Fig. ^ The con- 
tainer consists of an aluminium ring (diameter D —18 
cm) seperating two parallel glass windows by a gap of 

10 mm. It is filled by two unmixible liquids: a silicone 

011 (SOIL, Dow Corning, viscosity ij = 20 mPas, density 
p = 949 kg/m'^) and an aqueous solution of sugar and 
NiS04(WSS, Tj = 7.2 mPas, p = 1185 kg/m^). The liq- 
uid liquid interfacial tension has been determined with 
the pending droplet method to 35 ± 2 dynes/cm. The 
ratio of the filling heights SOIL/WSS was 8.4/1.6. The 
choice of heights was made in order to obtain a variety of 
different patterns, including a transition from squares to 
hexagons [3 ■ The sugar concentration has been adapted 
to match the refractive index to that of the covering sil- 
icone oil {n ~ 1.405) to a precision of 5 x lO^"'. The 
Ni^+-ions produce a broad absorbtion band in the spec- 
tral region 600-800 nm and provide high contrast pat- 
terns projected onto the diffusive screen. The container 
is illuminated from below with parallel light, and a band 
pass filter in front of the camera was used to detect only 
wavelengths A = (655±5)nm. By varying the intensity of 
the lamp the fiat interface has been set to a level of about 
50% of the maximum optical transmission. At a NiS04 
concentration of 17% by weight the contrast between the 
light intensity passing through crests and valleys of the 
wave pattern was optimum. The associated coefficient of 
optical absorption was measured as a = 5.2 ± O.lcm"^. 
In order to avoid uncontrolled changes of the viscosity, 
density and interfacial tension of both liquids all the 
measurements were performed at a constant temperature 
(23±0.1°C). The Faraday waves were excited by an elec- 
tromagnetic shaker vibrating vertically with an acceler- 
ation in the form a{t) ~ aQCOs(flt). The driving signal 
came from a computer via a D/A-converter and the ac- 
celeration has been measured by piezoelectric sensor. A 
self developed closed-loop algorithm was used to suppress 
higher harmonics n(Qt) in the driving signal to guaran- 
tee a purely harmonic driving. Faraday patterns were 
recorded in the following way: a high speed (250 Hz) 8- 
bit CCD camera was mounted some distance above the 
diffusive screen. Pictures were taken synchronous to the 
external driving. For a certain instant to the surface el- 
evation of the Faraday patterns h{x,y,to) is given by: 



h{x,y,to) = -In— — 

a Ip[x,y,to) 



(1) 



where Ir{x,y) and Ip{x,y,to) are 2D intensity distri- 
butions captured by the camera for the reference pic- 
ture (fiat interface, Uo = 0) and for the Faraday pat- 
tern (flo / 0), respectively. Finally, the surface elevation 
function h{x, y, t) are Fourier transformed and the time 
evolution of the Fourier amplitudes and phases of spa- 
tial modes is extracted. The use of a high speed camera 
compared to the earlier measurements by some of the co- 
authors allows for a better temporal resolution and the 
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FIG. 1: Experimental setup: L - halogen lamp, Ln - lens, C - 
container filled by two liquids: SOIL and WSS with the same 
refractive indices, DS - diffusive screen, IF - interference filter, 
CCD - high speed CCD camera. 



method is not sensitive to distortions (defects) on time 
scales of several periods. The logarithm of the intensity 
profile renders the dynamic range nonlinear, and with 
an 8-bit dynamical range the resolution is approximately 
1% (2%) at small (high) surface elevations, relative to 
the maximal surface heights. The validity of the method 
has been also checked with fiat layers of colored liquids of 
different thicknesses and the same accuracy was found. 
However, one should note the the Fourier transformation 
integrates over many pixels and a significant better res- 
olution is to expect. 



3. The linear regime 

The experiments have been performed by quasistati- 
cally ramping the driving amplitude for the frequencies 
/ = fj/27r = 12, 16, 20, 29 and 57iJz from shghtly below 
the critical acceleration Oc (e = (a — ac)/ac = —0.02) up 
to just below the acceleration where the interface disin- 
tegrates and droplets form. In the same form a ramp was 
driven down to check for hysteretic effects, of which none 
were found. For each amplitude step a series of pictures 
where taken and Fourier transformed. Typically the pat- 
tern occurs in the center region of the container first but 
evolves in a range of Ae = 0.02. From the Fourier trans- 
formation of the pictures at —0.02 < e < 0.1 (Fig. O 
the critical acceleration Qc and the critical wavenumber 
kc has been determined (Fig. [21). The experimental data 
can be compared with the results from the theoretical 
linear stability analysis that has been performed using 
the algorithm proposed by Kumar and TuckermanQ. 
The agreement between theory and experiment is very 
good, s imilar to former studies at the liquid-air interface 
haM. But with our new technique we are now able to 
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Equation n then reads 




N 
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driving frequency Ql2n (Hz) 
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FIG. 2: Critical acceleration ac and critical wave number kc 
for different driving frequencies fl. The symbols mark experi- 
mental data, the lines the theoretical linear stability analysis. 
The size of the symbols coincide with the size of the error 
bars. 



verify more details of the predictions of the linear theory, 
e.g. the temporal spectra at onset of the instability. It 
is a particular feature of the Faraday- Experiment, that 
at onset only one wave number kc becomes unstable, but 
the temporal spectrum already contains multiples of the 
fundamental oscillation frequency uj at onset. We arc in 
the regime of subharmonic response and the fundamen- 
tal oscillation frequency at onset is always ui — n/2 but 
the spectrum also contains {n+ 1/2)17 frequency compo- 
nents. More precisely we can write the surface deforma- 
tion h{r,t) as 



1 



N 



h{v,t) = - (A.e"^'"- + c.c.) 

i=l 



+ 00 

E 

n— — oo 



\(n+l/2)nt 



(2) 



Here r = (a;, y) are the horizontal coordinates. The set of 
complex Fourier coefficients are the components of the 
eigenvector related to the linear stability problem and de- 
termine the subharmonic time dependence. The spatial 
modes are characterized by the wave vectors , each car- 
rying an individual complex amplitude Ai. These quanti- 
ties are determined by the nonlinearities of the problem. 
In principle the k.; can have any length and orientation 
but at onset the relation |k,;| — kc holds. The number 
N of participating modes determines the degree of rota- 
tional symmetry of the pattern: = 1 corresponds to 
lines, = 2 to squares, N = 3 to hexagons or triangles, 
etc. It can be shown |23 that the C,n and C_„ are cou- 
pled in a way that = so that heterodyning of right 
and left travelling waves always result in standing waves. 



/i(r,i) =^(|A,|cosk, ■T^4>i) 

i=l 

x"^{\Cn\cos{n + l/2)nt + i:n) (3) 

n=0 

The complex eigenvectors Cn can be calculated modulo a 
constant factor and the ratio of the amplitudes |Cn l/ICn-i-! I 
as well as the temporal phases ipn can be compared with 
experimental data. They are obtained in the following 
way: For each step in the driving amplitude a series of 
snapshots of the surface state (Fig. 13)) is taken. The 
primary pattern consists of squares and their formation is 
governed by the nonlinearities of the problem and one of 
our goals is to identify how far from onset the predictions 
from the linear theory hold. 

An analysis of the Fourier transformation of the pic- 
tures yield amplitudes A(i)(k(zj)) that are shown in Fig. 
El For the wave vectors k(ij) the nomenclature from crys- 
tallography is used, e.g. k(10) and k(01) are the vectors 
that generate the simple unit cell of the square pattern 
(Fig. O. The temporal evolution of the amplitude of 
one of the critical modes A(i)k(10) with |k(10)| = kc 
(Fig. is then again Fourier transformed and a typi- 
cal spectrum is shown in Fig. [7^. These data are taken 
for all driving strengths e (Fig. |S^) and we always find 
the same values for A(t)k(10) and A(t)k(01) within the 
experimental resolution. In agreement with former inves- 
tigations .15,] in a system with a larger aspect ratio (con- 
tainer size to wave length) our study reveals also that the 
fundamental spatial mode |k(10)| = kc for all e. We can 
now extract the ratio of yl(fi3/2, k(10))/A(f]/2, k(10)) 
that is shown in the inset of Fig. |3| The contribu- 
tion of higher harmonics is in the order of 5 to 10% 
and increases slightly with the driving strength. The 
agreement between the experimental data and the lin- 
ear theory is again very good up to driving strength 
e > 0.5, especially for lower driving frequencies. It is 
very surprising that the agreement even holds up to sec- 
ondary patterned surface states, where a transition from 
a square to a hexagonal state has been taken place and, 
as we show later, strong nonlinear contributions partic- 
ipate in the dynamics of the system. The experimental 
data show also that at driving strength as low as e = 0.02 
the surface state consists of no measurable higher spatial 
Fourier modes (see Fig. ISJ but of higher temporal har- 
monics in perfect agreement with the linear theory. This 
allows an extrapolation of A{fl3/2, k(10))/A(fi/2, k(10)) 
to the neutral situation e = for all driving frequencies 
(Fig. 13) . The frequency ratio decreases first with increas- 
ing frequency and has a minimum at f2/(27r) « AQHz. 
This characteristic shape reflects the amount of damp- 
ing present in the system. At low driving frequencies 
the ration between fill hight and wave number is small. 
In this regime damping from the bottom, that increases 
with decreasing frequency, is most significant. For larger 
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FIG. 3; The ratio of the amplitudes A{m/2) / A{Q./2) of the 
k(10) mode at e = for different driving frequencies. The 
values are extrapolated from measurements at e > shown 
in the inset: the amplitude ratios at Q./2n — 12Hz and 29Hz 
as a function of the driving strength e. 



FIG. 5: Snapshots of the surface state and the power spectra 
at Q/2tv = 12Hz and e =0.17 (ao =30.0 m/s^) for two differ- 
ent temporal phases a) at maximum and b) minimum surface 
elevation as indicated in Fig. El 
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FIG. 4: The temporal phases ip {Q./2) and '4){ZQ./2) of the 
k(10) mode for two driving frequencies versus driving strength 
e. The symbols mark experimental, the lines theoretical data. 
Squares and broken line: Q,/2n — 57Hz. Circles and dotted 
line: Q/2tt = 12Hz; in the range 0.2 < £ < 0.28 a transition 
from squares to hexagons takes place and in this disordered 
state an extraction of phases is not possible. 



driving frequencies the damping from the bulk of the 
liquid (a function increasing with the frequency) is the 
strongest contribution. This behavior is also reflected in 
the critical accelerations (compare with Fig. |2J|. The 
ratio A(f^5/2,k(10))/yl(f]/2,k(10)) has been evaluated 
too, but the experimental resolution is not sufficient here 
for a conclusive comparison between theory and experi- 
ment. 

In the same way the temporal phases ■(/'„ can be ex- 
tracted from the Fourier spectrum and once more a good 
agreement between the theoretical predictions and ex- 
perimental data is obtained, at least for the fundamental 
ri/2 component. For the 3il/2 component the scatter of 
the experimental data is very large and we find significant 
differences between experiment and theory. But besides 
the large scatter we observe a pronounced nonmonotonic 
behavior of the ip3/2n component at 57 Hz and this part 
of the spectrum seems to be governed by nonlinear inter- 
actions. 



4. The nonlinear surface state at = 12Hz 

A. The square state 

The primary pattern near onset (0 < e < 0.28) con- 
sists of squares, shown in Fig. [S] Their formation is 
determined by the minimum of the Lyaponov functional 
of the according amplitude equation of the critical modes 
|2ll | and a quantitative theoretical prediction of the ex- 
pected pattern can be given by inspection of the cubic 
coupling coefficient Q. To our knowledge, for a two liq- 
uid system there has not yet been an attempt to calculate 
this coefficient, but squares are a common pattern in free 
surface experiments with low viscous liquids. The am- 
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FIG. 6: The absolute amplitudes A of diflterent spatial modes 
and the driving signal a{t) in the square state at Q/2tt — 
l2Hz and e =0.17 (ao =30.0 m/s^). Please note that un- 
like in Ref. not the square root of the power spectra but 
the amplitude A of a deformation h{r,t) = ylcos(k(y) ■ r) is 
shown. 



plitude equations follow from a solvability condition of a 
weakly nonlinear analysis of the underlying constitutive 
equation, and its principal form is determined by the 
symmetries of the system. For the subharmonic response 
one can write 



N 



,)|A(k,)|M(k,), (4) 



with r the linear relaxation time and T{9ij) the cubic 
coupling coefhcient that depends on the angle 9ij between 
the modes k^ and k^ with |kij [ = kc- The amplitudes 
A{ki) are modulated with a subharmonic (n+ 1/2)17 time 
spectrum given by the Qn from the linear eigenvectors. 
Equation 2] predicts a pitchfork bifurcation and in order 
to study this scenario one has to extract the different 
temporal Fourier modes of the measured A(ki, t) (Fig. O 
first. The result is shown in Figs. [TlandlSJ^. As long as 
the pattern consists of squares there are no harmonic time 
dependencies in the basic spatial modes to observe, but 
a continued growth of and 317/2 contributions. The 
512/2 contribution is very weak and only slightly larger 
than the noise. The square of the sum of the amplitudes 
As = A(r2/2) + A(317/2) yields a straight fine if plotted 
versus the driving strength e (Fig. as one would 

expect for the case of a pitchfork bifurcation. From the 
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FIG. 7: The temporal spectra of the amplitudes A(}i.(ij, uj)) 
of different spatial modes at Q./2ii = 12Hz and a driving 
strength: a) £=0.18 in the square state and b) e =0.38 in the 
hexagonal state. 



slope we can extract the cubic coupling coefficient Ag = 
e/r(90°), and we find r(90°) = 0.179mm-^. 

Now we can inspect the next higher harmonic spatial 
modes 74(k(ll)) and A(k{20)). Their temporal evolution 
is shown in Fig. El Both modes are a result of an inter- 
action of two fundamental modes, k(ll) — k(10) -I- k(01) 
and k(20) = k(10) + k(10). Quadratic coupling does not 
appear in the amplitude equations, but they are a nat- 
ural consequence of nonlinear spatial wave interaction 
and it is no surprise that we find that they obey har- 
monic oscillations, shown in Fig. [7^. The striking result 
of our analysis is rather the constant offset that we find 
in the A(k(ll)) and A(k(20)) spectrum (Fig. EJ. This 
means that in addition to the oscillatory part, the interfa- 
cial profile is also composed of contributions of constant 
deformations of the form h{r,t) — |Ai|coski-r. This 
might surprisingly first, but please note that this does 
not violate the mass conservation. Actually, it is a sim- 
ple consequence of the quadratic coupling of a real stand- 
ing surface wave oscillation, 3fi(e'''* '"e'''* ''e*^/^e~'^/^) — 
cos2ki ■ r(l -t- cosil) (compare also Fig. 114b'). 



This quadratic coupling scheme can be verified by plot- 
ting A(k(20)) and yl(k(ll)) versus the square of the am- 
plitude of the fundamental mode A(k(10)) or the product 
A(k(10)) X A(k(01)) respectively. The data can be per- 
fectly reproduced by a linear fit. From the slope one gets 
the strength of this nonlinear coupling and we do find 
the same values for all frequencies ll/27r — 16,20,29 Hz 
where squares are to be observed. Finally our Fourier 
analysis yields that the imaginary part of the coupling 
scheme obeys the same resonance conditions, and the 
spatial phase of the higher harmonic modes is given by 
0(k(2O)) = 20(k(lO)) and 0(k(lO)) -I- (/'(k(Ol)). 
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FIG. 8: The amplitudes A{n/2Q) of the a) k(10) and b) k(ll) 
mode at fl/2n = 12Hz as a function of the driving strength 
e. c): the square of the sum of the subharmonic components 
of the A{10) mode versus e. As expected for a forward bifur- 
cation the data can be hnearly fitted, at least up to driving 
strength of e ~ 0.1 
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FIG. 9: The amplitudes A(k(20)) and A(k(ll)) versus the 
square of the amplitude of the fundamental mode yl(k(10)) or 
the product of y4(k(10)) and yl(k(01)) respectively (f7/27r = 
12Hz). 



B. The hexagonal state 

In the range (0.20 < £ < 0.28) the pattern becomes 
disordered and transforms at higher driving strength to 
a hexagonal state (see Fig. I10|l that consists of three fun- 
damental spatial Fourier modes ki_2,3. But please note 
that for the construction of the crystallographic simple 



FIG. 10: Snapshots of the surface state and the Power spectra 
at fl/2-K = 12Hz and e =0.37 (ao =39.3 m/s^) for three 
different temporal phases, a) down hexagons, b) minimal 
surface elevation, c) up hexagons. See Refs. [3. ll2l| for further 
explanations on the switch from up to down hexagons in the 
Faraday- Experiment. 




FIG. 11: Vector diagram of the interacting modes for the 
hexagonal surface state. 



unit cell two vectors k(10) = ki and k(ll) = k2 are suf- 
ficient (k(il) + k(10) = k(01) = ka, as indicated in Fig. 

EH. 

An analysis of the temporal behavior of the amplitudes 
A{'k{ij)) of the spatial modes reveals, besides the the 
striking offset with the according constant spatial sinu- 
soidal surface deformation, both harmonic and subhar- 
monic time dependencies (see Fig. While harmonic 
(nfl) temporal contributions in the higher spatial har- 
monics k(20, 11) appear in a similar manner to that of the 
square pattern, the resonance between k(10) -I- k(Oi) = 
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lation of the 3D problem and quantitative tests for future 
work arc most important. In this sense we would like to 
encourage such attempts. But with our technique we 
are also able to verify known predictions from the lin- 
ear stability analysis and we find good agreement up to 
high driving strength of £ w 0.5. In the nonlinear state 
the most pronounced result is the identification of strong 
temporal constant sinusodial surface deformations in the 
spectrum. And with our possibility to access any Fourier 
component separately we can identify several resonance 
mechanisms, including an interesting case of a temporal 
resonance violation by use of spatial resonances. 
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FIG. 14: a,b): The absolute amplitudes A of different spa- 
tial modes in the line state at f2/27r =57Hz and e =0.11 
(ao =116.3 m/s^). In b) the temporal constant offset A{2k, 0) 
is indicated by the dotted line, c) The amplitude A{2k, 0) ver- 
sus A'^{k,Q/2) 
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